An experimental proposal for revealing contextuality in almost all qutrit states 
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Contextuality is a foundational phenomenon underlying key differences between quantum theory 
and classical realistic descriptions of the world. Here we propose an experimental test which is 
capable of revealing contextuality in all qutrit systems, except the completely mixed state, provided 
we choose the measurement basis appropriately. The experiment is based upon a recently derived 
result [Phys. Rev. A 86, 042125]. The 3-level system is furnished by the polarization and spatial 
degrees of freedom of a single photon, which encompass three orthogonal modes. Projective mea- 
surements along rays in the 3-dimensional Hilbert space are made by linear optical elements and 
detectors which are sensitive to single mode. We also discuss the impact of detector inefficiency and 
losses. 
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I. INTRODUCTION 



In classical physics, all observable quantities have an objective reality. The outcome of measuring a classical 
observable A can not be influenced by a property of the observer or by the choice to simultaneously co-measured A 
with a second observable B. This property is referred to as non-contextual realism. 

Distinctively, in quantum theories the measurement outcome for an observable A depends on what you choose to 
co-measure at the same time. There exist certain sets of observables where it is impossible to simultaneously assign 
pre-existing outcomes. In this scenario there is no joint probability distribution from which the measurement statistics 
for every observable can be recovered as marginals. 

Recently Klyachko, Can, Binicioglu and Shumovsky introduced a 5-ray (KCBS) inequality [T|. The KCBS inequality 
uses the maximum attainable level of correlation between certain co-measurable observables belonging to any non- 
contextual realistic theory. By construction the KCBS inequality is satisfied by all non-contextual hidden variable 
models. However Klyachko et. al. discovered violation of the inequality in a 3-level quantum system. It has been 
shown that the KCBS inequality is the simplest possible test for contextuality of quantum theories, in the sense that 
any inequality based on projective measurements along less than five different direction will be satisfied by all classical 
and quantum theories [2J. The KCBS inequality formed the basis of several recent experimental tests of quantum 
contextuality. Lapkiewicz et. al. have experimentally demonstrated violation of the KCBS inequality in an indivisible 
3-level quantum system (furnished by a single photon) [3] . This experiment built upon a robust and rapidly developing 
body of experimental continuality literature [4-9 . In our present work we draw on the tests outlined in these papers. 

From the perspective of an experimental test of non-contextual realism the KCBS inequality is a state dependent 
test; for a sufficiently pure state there exists a choice of five measurement directions which can be used to detect 
violation of the KCBS inequality. Hence any experimental implementation requires a rigorous state preparation 
phase. For 3-dimensional systems there also exists a completely state independent test involving thirteen projectors 
developed by Oh and Yu |10j and an almost state independent test involving nine projectors The latter is capable 
of revealing contextuality in every state, except the maximally mixed state, provided we choose the measurement basis 
appropriately. We describe this inequality in more detail in Section |IH] It has been shown that thirteen rays is the 
minimum number of distinct projector directions necessary for a state independent test of non-contextual realism [12] . 

The purpose of this paper is to propose an experimental test for the inequality from Ref. [llj . The experimental 
set up will be along the lines of Lapkiewicz et al.'s KCBS inequality experiment [3J. We will also discuss the minimum 
detector efficiency required to ensure the results arc independent of any auxiliary assumptions about photons in 
undetected events. 

In Section [H] we discuss the key differences between non-contextual hidden variable theories which obey classical 



realism and quantum mechanics. We introduce the inequality and contextual results from Ref. |llj in Section III In 
the following section we propose an experimental test of this inequality. We use Section [V] to analyze the effects of 
losses and detector inefficiency on this experiment. We conclude the paper with a brief summary of our results. 



II. CONTEXTUALITY 



Consider a theory with n observables, A\, . . .A n , where each observable, Aj, returns outcome a, with probability 
p(Ai — a,). This theory exhibits classical realism whenever there exists a joint probability distribution for the 
combined measurement outcomes p(A\ = oi, . . . , A n = a n ) such that the measurement statistics, p(Ai = ai), of each 
observable, Aj, can be obtained as a marginal of this distribution [13] . In this scenario there exists a non-contextual 
hidden variable model, conditioned on a set of parameters A. The joint probability distribution is reproduced through 



p{Ai = a%, . . . ,A n = a n ) = J dX p(X)p(A 1 = oi|A) . . .p{A n = a n \X), (1) 

where each strategy for assigning measurement outcomes, A, has weighted probability p(X) on the set A, and p(Aj = 
di\\) is the conditional probability of obtaining an outcome a, when measuring observable A4 (for the given non- 
contextual hidden variable A). 

Quantum theories do not necessarily exhibit non-contextual realism. When we measure an observable A\ we 
simultaneously collapse the state we are observing onto an eigenstate of A\. If A\ is independently simultaneously 
co-measurable with either member of a non-commuting pair [A% and A3), then choosing to measure A\ simultaneously 
with A2 (respectively A3) can erase information about the correlations between A\ and A3 (respectively A2). Because 
a joint measurement of A\ and A 2 can erase different information from a joint measurement of Ai and A 3 , the two 
joint measurements do not need to be unitarily equivalent. The choice to simultaneously co-measure A\ with either 
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Ai or A3 can influence the outcome of Ai. We say that these two non-commuting observables each furnish a context 
for A\ and that quantum mechanics is a contextual theory. 

Given a set of observables we can succinctly represent which subsets are compatible in a single graph. This graph 
has one vertex, Wj, for each observable A4. The edge set of the graph is generated by connecting any two vertices Vi 
and Vj associated with compatible observables, respectively Ai and Aj, by an undirected edge (vi,Vj). 



III. A SPECIFIC TESTS OF CONTEXTUALITY 



In this section we will explain a recently derived inequality from Ref. [11] , This inequality is based on nine 
projective measurements whose compatibility relations are provided by orthogonality and depicted in Figure [T] The 
nine vertices in Figure [I] can be identified with projective measurements II,; = along the following nine rays in 

a 3-dimensional Hilbert space: 

1) = (1,0,0), 

2) = (0,1,0), 

3) = (0,0,1), 

4) = 1/72(0, 1,-1), 

5) = 1/73(1,0,-72), 

6) = 1/73(1,72,0), 

7) = 1/2(72,1,1), 

8) = 1/2(72, -1,-1), 

9) = 1/2(72,-1,1). (2) 

In a non-contextual theory, for any complete orthogonal basia^of rays it is possible to assign values of 1 or 

to the corresponding projectors II ni , in such a way that Yi u J\ Ui = 1 (completeness of the basis) and if a projector, 
H Ui , appears in two different complete basis then it is assigned the same consistent value in each basis (non-contextual 
condition) [T4"] . 

It follows that for any non-contextual realistic theory, where measurement statistics can be encapsulated in a hidden 
variable strategy p(X) as in ([T]), these nine-projectors satisfy: 



(3) 



i=i 



We can arrive at the contextuality inequality in Ref. by replacing each of the projectors in Equation ([3]) by a 
dichotomous observable, 



Ai = l — 2tti, 

which assign values of +1 and —1 to measurement outcomes. In this notation Equation d3j) 
about the correlations of observables A\,...,A$ (in relation to the edge set E(G) of Figure [ll 



(4) 



becomes a statement 



(5) 



Given a quantum state p, which is not the maximally mixed state, it is possible to choose a basis for the 3- 
dimcnsional Hilbert space so that the inequality (|3| will be violated. More explicitly, the state p violates Inequality 
([3]) when you make projective measurements along the directions in Equation ^ relative to this basis jllj . Optimal 
violation will be obtained when we choose to measure with respect to the eigenbasis of p. To complement this analysis 
we have run a numerical simulation, using QI Mathematica package |15j . which tested the inequality ^ for a prefixed 
measurement basis and a randomly generated sample of ten million density matrices weighted by the Hilbert-Schmidt 
measure. We found 49.98% of quantum states p violated 



For the d-dimensional Hilbert space. The pertinent case being d = 3. 



FIG. 1. The graph, G, of the orthogonality relations for the nine projectors in the inequality from Ref. 



For the purpose of the following experimental design it is pertinent that the projectors |1), |2) and |3) form a 
complete basis for the 3-dimensional Hilbert space. Hence for these projectors there is a completeness relation: 

(AxA 2 ) + (A 2 A 3 ) + {A Z A X ) = -1. (6) 

A similar relation holds for the complete basis |7), |4) and |8). 

IV. EXPERIMENTAL DESIGN 

Lapkiewicz et al. have experimentally demonstrated violation of the KCBS inequality for the indivisible 3-level 
quantum system furnished by a single photon [3j. A photonic qutrit uses two distinct spatial paths for the photon 
and within one of the spatial paths two optical modes propagate as distinct polarization modes. Here we propose 
a straightforward adaptation of Lapkiewicz et al.'s experiment to the inequality ([5]). The aim of our proposed 
experiment is to directly measure violation of the inequality ^ by a 3-level quantum system. This will involve taking 
simultaneous measurements of pairs of compatible observables Ai and Aj when the corresponding vertices Vi and Vj 
are connected by an edge in Figure [T] An implicit caveat on any experimental implementation of this proposal is 
that all experimental runs must be independent and the results after a statistically significant number of runs must 
represent an unbiased sampling from the joint probability distribution ([!]) ; in the case where the latter is tailored to 
the nine projective measurement constructed from ([2]). Later we will discuss how losses and detector inefficiencies may 
allow a joint probability distribution conditioned on a hidden variable to violate the inequality ([5]). In both scenarios 
we use a single photon heralded source which will alow us to monitor photon losses. Previous experiments, 0[5], 
produce pairs of polarization entangled photons in the singlet state via parametric down conversion and subsequently 
rerouted one photon to an auxiliary detector 'detector 0'. Post selecting on whether this auxiliary detector clicks 
creates a heralded single photon, see Figure [3] Any desired initial state can be synthesized using polarizing beam 
splitters to combine the two spatial modes and half wave plates to transform the polarization components inside a 
single spatial mode. 

We introduce a single detector for each mode (we will use the label: detector i, for i = 1,2 or 3). Two scenarios are 
considered: in the first scenario there are no losses and detector i will click whenever a photon is in the corresponding 
mode In the second scenario we will allow for inefficient detectors. We revisit this scenario later. It is convenient 
to identify the three orthogonal modes with rays |1), |2) and |3) from Equation ([2]). If detector i clicks then we record 
H; = 1 and using Equation ^ we assign a value Ai = — 1 to the corresponding dichotomous observable. 

The correlations of two dichotomous observables are given by: 

(AiAj) = P(Ai = A, = 1) + P(Ai = A, = -1) - P(Ai = 1, Aj = -1) - P(A, = -1, Aj = 1) 

= 1 - 2P(A, = 1, Aj = -1) - 2P{Ai = -1, Aj = 1). (7) 

The second line follows because (with unit probability) either both Ai and Aj are equal to 1, alternatively both are 
equal to —1, or precisely one of them is equal to —1. We only count results from experimental runs where a single 
detector clicks coincidentally with detector 0. This implies that in all events contributing to Equation Q the sum of 
the measurement outcomes for the dichotomous observables (during a single run) is equal to —1 (in the first instance 
Ai + A 2 + A 3 — — 1). Furthermore when Ai and Aj correspond to orthogonal modes we discount experimental runs 
where Ai = Aj = —1 causing P(Ai = Aj = —1) = 0. 
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Repeated runs of this experimental configuration can be used to measure the correlations (A1A2), (A1A3) and 
(A2A3). A full experiment capable of measuring all the correlation in the inequality ([5| can be obtained by adding 
standard optical elements; specifically half wave plates and polarizing beam splitters. If we pass the input modes 
|1), |2) and |3) through a series of half wave plates and polarizing beam splitters then we can output the orthogonal 
modes for any three rays \j) and \k) which are a linear combination of |1), |2) and |3), and form a basis for the 
3-dimensional Hilbert space. At every stage of the experiment each detector is aligned to detect photons in a single 
mode. It is assumed that this detector will click if the photon is in the corresponding mode. 

A schematic of the sequence of optical elements and measurements needed to obtain all the correlations in the 
inequality |5| is depicted in Figure [2j The physical implementation of this experiment will require exact details 
regarding the sequence of half wave plates and polarizing beam splitters needed to appropriately mix the two different 
polarization modes within a single spatial mode and subsequently mix the spatial modes. We leave the nuances of the 
setup to the discretion of an experimental group. In Figure [3] we give a proof of principle that the schematic given in 
Figure [2] is physically realizable using an adaptation of Lapkiewicz et al.'s design [3J. For more detail see Appendix 
[X] which contains Figure [3] and lists the transformations enacted by the half wave plates. 




FIG. 2. Each subfigure represents an experimental configuration which measures (AiAj) for the specific Ai and Aj indicated 
on the righthand side of the subfigure. Orthogonal modes are represented by horizontal lines. Sequences of half wave plates 
and polarizing beam splitters are labeled T\, . . . , Ts and detector positions are indicated by the corresponding dichotomous 
observable. The expectation value (Ag) can be obtained from the data collected during the experimental runs depicted in 
subfigures 6) or 7). The vertical line through the detectors in subfigures 1) and 3) indicate it is not necessary to record data 
during these stages. This follows from Equation (IM. The precise sequence of optical elements Xi, . . . , Tg are recorded in Figure 
[3]which should be read in conjunction Table \K\ in Appendix \K\ In this schematic when two consecutive stages of the experiment 
involve measuring the same mode, \i), and the beam line of mode \i) is unobstructed between the points at which these two 
measurements are made, we assume the physical implementation will measure the (same) mode \i) in two different contexts. 



The difference between A2 in subfigure 1) and A' 2 in subfigure 5) is explained in detail in Section IV 



It is implicitly assumed that at any stage of the experiment; if a given mode's beam line is not interrupted by an 
optical element then there is perfect transmission of any photon (amplitude component) in this mode. This occurs 
regardless of whether the other two modes pass through a half wave plate or are obstructed. 

Notice that according to Figure [2j we first measure the mutual correlations of the dichotomous observables, 
Ai,A 2 ,A 3 , corresponding to orthogonal modes |1), |2) and |3). Subsequently in subfigure 5) we need to measure 
(A2A5), however during the intermediary stages the beam line of mode |2) has been interrupted. This introduces 
a source of error. We must recreate the mode |2) before measuring (A2A5). We can not be certain that in both 
instances we are measuring the precisely the same orthogonal mode |2). To accommodate this we distinguish between 
the observable A2 in subfigure 1) and the observable A' 2 measured in subfigure 5) of Figure [2j We correspondingly 
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modify the inequality to 

[(AiAj)(l - 6 l5 6 :j2 - S it6 S j8 ~ 6 i6 Sj 3 )] + (A 5 A' 2 ) + {A & A'%) + (A 6 A 3 ) + (A 9 ) >-l-e a -e b - e c , (8) 

(i,j)eE(G) 

by including the error terms e a = 1 — (A 2 A' 2 ) , e b = 1 — (A S A' S ), e c = 1 — (A 3 A' 3 ). 

In principle the method for obtaining each error terms is identical and is outlined in Ref. [3j. Therefore without 
loss of generality we explicitly focus on measuring the error term e Q = (A 2 A' 2 ) . We can obtain e a directly by using 
Equation ^ to determine the value of (A 2 A' 2 ). It is not possible to measure modes |2) and |2') concurrently. However 
we conclude from Equation Q that: 

e a = 2P(A 2 = l,A' 2 = -1) + 2P(A 2 = -1, A 2 = 1). (9) 

To find the above probabilities we set up the experiment with the intention of measuring A' 2 but place a carcj^] at 
the point where mode |2) is first created. This card should effectively block photons belonging to second orthogonal 
mode in sub figure 1) of Figure [2j This will exclude all events where A 2 = — 1 . If the detector measuring A' 2 clicks 
then we record an event where A 2 = 1 and A' 2 = —1. This will allow us to find P(A 2 = 1,A 2 = —1). We can chose 
to unblock mode |2) and instead block the other two modes: |1) and |3) in Figure ^1 In this scenario we exclude all 
events where A 2 = 1. When the detector measuring A' 2 does not click we have an event where A 2 = — 1 and A' 2 = 1. 
This will allow us to find P(A 2 = -1, A 2 = 1). 

Finally we highlight that Equation ([6| implies we do not need to collect data during the experimental stages 
depicted in subfigures 1) and 3) of Figure [2} We have drawn a vertical line through the detectors whenever it is not 
necessary to collect data. We include these stages in the schematic because it is important to comprehensively set 
up all the measurement configurations. We need to ensure the orthogonal modes have the same relative orientations 
(and orthogonality relations) as the projectors in Equation This makes it essential to setup all stages of the 
experiment. 

Errors will arise from imperfect alignment of half wave plates. This will cause some of the modes to be be misaligned 
with the rays in Equation ([2]). In this case we will be testing the inequality ^ using projective measurement along 
directions which do not perfectly match the list in Equation ([2]). The nine direction in Equation ([2]) corresponded to 
a theoretically optimal set of projective measurements for testing the inequality |5]). The errors introduced by half 
wave plate misalignment will make it more difficult to experimentally violate this inequality. However they will not 
invalidate any experimental results which directly demonstrate violation. 

We now come back to the second scenario: in any realistic experiment there will be runs where the photon is lost 
and none of the detectors click. These losses are caused by reflections off half wave plates and detector inefficiencies. 
To ensure the results will not be biased by losses we only consider data from experimental runs where a single detector 
clicks. This approach requires the subsidiary assumption that for all runs the probability of loosing a photon is not 
correlated with the measurement outcome. In the next section we consider how efficient the experiment must be in 
order to eliminate this assumption. 



V. DETECTION LOOP HOLE 



Without loss of generality we can assume all photons are lost at the detector. We attribute an efficiency r] to each of 
the three detectors. We want to make our results independent of the distribution of photon modes in undetected runs. 
To do this we include the runs where none of the detectors click. The combined total number of runs is constituted 
by runs where either a single detector clicks or none of the detectors click. 

In the ideal scenario where there are no losses, p(Ai = ai) is the probability of getting outcome ai when you measure 
observable Ai. When we include losses we use p'(A; L = a,;) = i]p(Ai — a») for the proportion of the combined total 
number of runs where we measure Ai and get outcome a^. 

To uniquely specify the results of a run we must specify the outcomes ai,a.j and a k of all three dichotomous 
observables (respectively Ai, Aj and A k ). 

When estimating the correlations in Equation ([T]) we now use the renormalized probability: 



p'(A i = a i ,A j = a j ,A k = a k ) 
p(Ai = at, Aj = aj ,A k = a k ) = ^-^—^—^—^ , (10) 



2 A practical implementation of this experiment uses only two spatial modes. At each stage of the experiment two of the levels (in the 
three level system) will have the same spatial mode index (but opposite polarization indices). Therefore the practical implementation 
of this scheme may require the card to be opaque only to photons in the polarization mode associated with |2). 
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which is the ratio of the number of runs with the given outcome (Aj = di,Aj = aj,A k = dk) to the number of 
runs where a single detector clicks. This eliminates the dependence on no detection events (runs with outcome 
At = Aj = A k = 1). 

In hidden variable theories we can use ([T]) to express p(A^ — a>i,Aj — aj,A k = a,k) as a function of the detector 
efficiency r\ and probabilities in the ideal scenario where there are no losses: 



p(Ai = di,A 3 = a,j,A k = a k ) = - - — - — J dX p(X) ■qp(A l = a^X) rjp(Aj = a 3 \X) i]p(A k = a k \X). (11) 

This allows us to evaluate the correlations in Equation Q for a hidden variable strategy p(X) conditioned on A. We 
find the inequality ^ is modified in the presence of detector inefficiency to 

£ (AiAj) + (A 9 ) > -4/ry 2 . (12) 

(:hj)eE(G) 



We need r\ > 0.82 for the right hand side of Equation (12) to be below the maximum level of quantum violation 



for this inequality [TB]. Hence when the detector efficiency is below 82% it is impossible to experimentally detect 
violation of the inequality ^ unless subsidiary assumptions are made about the distribution of putative outcomes 
(which detector failed to click) for undetected events. 



VI. CONCLUSION 



In conclusion we propose an experimental test of the inequality from Ref. [11] . An experimental implementation 
would test one of the foundational principles which distinguish quantum mechanics from classical realistic descriptions 
of the world. Furthermore if this test reveals contextual behavior in an indivisible spin one system then the behavior 
can not stem from entanglement. This result complements work done on the KCBS inequality [2 [3] as well as the 
thirteen projector state independent inequality from Ref. [5l 110] . We present numerical studies which indicate the 
inequality ([5| is violated by 49.98% of states in a 3-level system. This makes an experimental test less sensitive to 
state preparation. 
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Appendix A: Half wave plate orientations 



This appendix contains details of a possible implementation of Schematic [2j 

The table lists the transformations that half wave plates WPl — WP6 enact on the two polarization modes. 
Each angle corresponds to the angle of rotation in a 2 x 2 rotation matrix acting on a pair of polarization modes 
within a single spatial path. We have chosen to give the angle characterizing the transformation rather than the 
precise configuration of the wave plate's optical axis because the latter will be contingent on the exact details of 
the experimental implementation and may be obsolete if even a small change is made. The '— ' indicates that the 
corresponding wave plate is unnecessary for measuring the corresponding correlation. 



a) Preparation 



LD ppKTP F PBS 

—HI ^ 






b) Measurement 



WP1 




FIG. 3. Here we give a possible experimental implementation of Schematic[2]in the spirit of Lapkiewicz et. al. [3]. Preparation: 
a single photon heralded source is created through parametric down conversion; a process which outputs pairs of polarization 
entangled photons in the singlet state. Subsequently a polarizing beam splitter is used to redirect one member of the pair 
to a heralding detector: 'detector 0' which is denoted by '0' in the figure. Post selecting on events where detector clicks 
clearly identifies runs where the photon is lost (detector is the only detector to click) as distinct from runs where the photon 
is not lost (one of detectors 1-3 click coincidentally with detector 0). We take inspiration from Ref. [3] which used a 405 
nm Laser diode (LD) in combination with a nonlinear-periodically-poled-Potassium-Titanyl-Phosphate-crystal (ppKTP) and 
suitable filters (F) to produce pairs of polarization entangled photons. For a comprehensive explanation of this process see 
Ref. j3]. Any desired initial state can now be synthesized using additional half wave plates and polarizing beam splitters. 
Measurement: there are always two concurrent modes with the same spatial index and opposite (orthogonal) polarization 
indices. Half wave plates WP1 — WPG manipulate these two modes. Each wave plate implements a 2 x 2 rotation on the 
two dimensional subspace spanned by the pair of polarization modes within a single spatial path. The angles characterizing 
these rotations are given in Table [A] Modes with different spatial indices can be recombined using polarizing beam splitters 
which are represented by grey rectangles. We use radians to denote the configuration of an optical element (alignment of the 
optical axis) which implements the identity transformation. The remaining light-blue wave plates have been inserted to balance 
the path lengths of the two spatial modes; these are all set to n/4 (this alignment balances a ir phase shift between the two 
spatial modes). When the run requires a subset of the optical elements, the unnecessary elements can be aligned to identity 
transformations. Note that in Figure [2] transformations T4 and T5 act consecutively on the same two-mode subspace. In the 
physical implementation T4 and T5 correspond to wave plate WPi set to two different orientations, see Table [X] Analogously 
T7 and T$ use two different settings of WPG. 
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